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Abstract 

We consider a hierarchy of the natural type Hamiltonian systems of n degrees 
of freedom with polynomial potentials separable in general ellipsoidal and general 
paraboloidal coordinates. We give a Lax representation in terms of 2 x 2 matrices 
for the whole hierarchy and construct the associated linear r-matrix algebra with 
the r-matrix dependent on the dynamical variables. A Yang-Baxter equation of 
dynamical type is proposed. Using the method of variable separation we provide 
the integration of the systems in classical mechanics conctructing the separation 
equations and, hence, the explicit form of action variables. The quantisation 
problem is discussed with the help of the separation variables. 



Key words: separable systems, linear r-matrix algebras, classical Yang-Baxter 
equation, Lax representation, Liouville integrable systems 

AMS classification: 70H20, 70H33, 58F07 



2 



1 Introduction 



The method of separation of variables in the Hamilton- Jacobi equation, 

dW 



H(p 1 , . . . ,p n , Xi, . . . ,X n ) = E, Pi 



, i = l,...n, 

OXi 



Xl) 



is one of the most powerful methods for the construction of action for the Liouville 
integrable systems of classical mechanics ||. We consider below the systems of the 
natural form described by the Hamiltonian 



H= - ^2pi+V(x h ...,x n ), Pi,XiGR. 
z i=i 



[1.2) 



The separation of variables means the solution of the partial differential equation ( |1 . 1|) 
for the action function W in the following additive form 

n 

W = y £W i { l n;H 1 , . . .,H n ), H n = H, 



where /ij will be called separation variables. Notice that the partial functions W{ depend 
only on their separation variables /ij, which define new set of variables instead of the 
old ones {xk}, and on the set of constants of motion, or integrals of motion, {Hk}. In 
the sequel we shall speak about coordinate separation where the separation variables 
{fii} are the functions of the coordinates {xk} only. (The general change of variables 
may include also corresponding momenta {pk})- 

For a free particle (V = 0), the complete classification of all orthogonal coordinate 
systems in which the Hamilton- Jacobi equation ( |1 . 1| ) admits the separation of variables 
is known, they are generalized n-dimensional ellipsoidal and paraboloidal coordinates 
||, |] (see also the references therein). It is also known that the Hamiltonian systems 
( |1.2| ) admitting an orthogonal coordinate separation with V ^ are separated only in 
the same coordinate systems. 

The modern approach to finite-dimensional integrable systems uses the language of 
the representations of r-matrix algebras jTDJ 15, |17| 18 . The classical method of separa- 
tion of variables can be formulated within this language dealing with the representations 
of linear and quadratic r-matrix algebras [11], 12, [17] 0. For the 2x2 L-operators 



the recipe consists of considerating the zeros of one of the off-diagonal elements as the 
separation variables (see also generalization of this approach to the higher dimensions 
of L-matrix (19[]). For V — in [11, |12| 2x2 L-operators were given satisfying the 



standard linear r-matrix algebra [10, 15 





( 1 










1 








1 





u 





1 










V o 








1 / 



;i.3) 



and the link with the separation of variables method was elucidated. Here in (|1.3| ) we 
use the familiar notations for the tensor products of L(u) and 2x2 unit matrix /, 
Lx(u) = L{u) ® J, L 2 (v) =1® L(v). 
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In the present article we construct 2x2 L-operators for the systems (|1.2|) being 
separated in the generalized ellipsoidal and paraboloidal coordinates. In the case when 
the degree N of the potential V is equal to 1 or 2 the associated linear r-matrix algebra 
appears to be the standard one ( |i.3p . In the case N > 2 the algebra is of the form 

{L^u), L 2 (v)} = [r{u - v), Ltiu) + L 2 {v)\ + [s(u, v), L^u) - L 2 {v)\, (1.4) 

with s(u, v) = ocn(u,v) cr_ <g> cr_, where a_ = o\ — ia 2 and Oi are the Pauli matrices, 
and cxn(u,v) is the function which equals 1 for N = 3 and depends on the dynamical 
variables {xk,Pk}t=i f° r N > 3. 

The study of completely integrable systems admitting a classical r-matrix Poisson 
structure with the r-matrix dependent on dynamical variables has attracted some at- 
tention flU H, [13]. It is remarkable that the celebrated Calogero-Moser system, whose 
complete integrability was shown a number of years ago (c.f. ||14j|), has been found only 
recently to possess a classical r-matrix of dynamical type |4| . 

Below we recap briefly how to get the 2x2 L-operators for the separable systems 
( |1.2| ) without potential V [11], [12| . Let us consider the direct sum of the Lie algebras 
each of rank 1: A — ©& = i sofc(2, 1). Generators Sk G R 3 , k = 1, . . . , n of the A algebra 
satisfy the following Poisson brackets: 

{sh s li} = s km£ijigiis{, = diag(l,-l,-l). (1.5) 

Everywhere below, we will imply the g metric to calculate the norm and scalar product 
of the vectors sf. 

s? =(4Si) = (^) a -(«?) a -(«?) 2 , 
{si, = s i s k — s 2 s| — s 3 s 3 . 

Let us fix the values of the Casimir elements of the A algebra: sf = c 2 , then variables 
Si lie on the direct product of n hyperboloids in R 3 . Let q G R and choose the upper 
sheets of these double-sheeted hyperboloids. Denote the obtained manifold as K+. We 
will denote by hyperbolic Gaudin magnet [0 integrable Hamiltonian system on 
given by n integrals of motion Hi which are in involution with respect to the bracket 

(ED, 

tf. = 2 £^3 {Hi,H k }=0, e^e k if i^k. (1.6) 

fe=l e * ek 

To be more exact one has to call this model an n-site so(2,l)-XXX Gaudin magnet. 
Notice that all the Hi are quadratic functions on generators of the A algebra and the 
following equalities are valid 



£i?i = o, = J 2 -]T 

i=l i=l i=l 



where the new variable J = Yh=i Si is introduced which is the total sum of the hy- 
perbolic momenta Si. The components of the vector J obey so(2,l) Lie algebra with 



respect to the bracket ( |1.5| ) and are in involution with all the Hi. The complete set 
of involutive integrals of motion is provided by the following choice: Hi, J 2 and, for 
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example, (J 3 ) 2 . The integrals ( |1.6| ) are generated by the 2x2 L-operator (as well as 
the additional integrals J ) 



L{u) = £ 



j=l u- e] \-{s) + s)) -is) )> 



det L{u) = - V —1- + —- 1—- 9 , (1.7) 




satisfying the standard linear r-matrix algebra ( |1.3| ). Let q = 0, i = 1, . . . ,n, then 
the hyperboloids sf = cf turn into cones. The manifold K£ admits in this case the 
following parameterization (pi, X{ G R): 

9 9 9 9 

1 _ Pi ~r x i 2 _ x i 3 _ A x » /-i q\ 

where the variables pi and x« are canonically conjugated. Using the isomorphism ( |1.8| ) 
the complete classification of the separable orthogonal coordinate systems was provided 
in [1TTL ITJ] by means of the corresponding L-operators satisfying the standard linear r- 



matrix algebra ( |1.3| ). The starting point for our investigation are these L-operators 
written for the cases of free motion on a sphere and in the Euclidean space. 

The paper is organized as follows. In Section 2 we describe the classical Poisson 
structure associated with the hierarchy of natural type Hamiltonians separable in the 
three coordinate systems — spherical (for motion on a sphere), and general ellipsoidal 
and paraboloidal (for n-dimensional Euclidean motion) coordinates. This structure is 
given in terms of the linear r-matrix formalism, providing new example of the dynamical 
dependence of the r-matrices. We also introduce analogue of Yang-Baxter equation for 
our dynamical r-matrices. In the Section 3 we derive the Lax representation for all 
the hierarchy as a consequence of the r-matrix representation given in the Section 2. 
Section 3 deals also with the variable separation material. The question of quantisation 
of the considered systems is briefly discussed. 

2 Classical Poisson structure 

Let us consider the following ansatz for the 2x2 L-operator 
where 

n 2 

B{u) = e-J2^^, e = 0,l, or 4(u-x n+1 + B), (2.2) 

A{u) = e + lt^, (2.3) 

2 u — ti 
t=i 1 

n 2 N 

Cn(u) = J2^^-V N (u), V N (u) = Y.V k u N ~ k . (2.4) 



i=l U e i k=0 
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Here the canonically conjugated variables ({pi,Xj} = Vk are indetermi- 

nate functions of the x- variables; B and e« are non-coincident real constants. Note that 
dot over e means differentiation by time, and for natural Hamiltonian ( |1.2|) one has 

%n+l = Pn+1- 



orm 



Theorem 1 Let the curve det(L(u) — XI) = for the L-operator ( \2. 1\ ) have the f( 
A 2 - A{uf - B(u)C N (u) = \ 2 + eu N -Y = 0, 

~i U - Ci 

for e = 0,l, and (2.5) 

n tt 

A 2 -A(uf -B(u)C N (u) = \ 2 + l6u N - 2 (u + B) 2 + 8H -Y — = 0, 

for e = 4(w - x n+ i + B) , (2.6) 

with some integrals of motion Hi and Hi, H in the case of ( $.(\ ). Then the following 
recurrence relations for Vk are valid 

n k— 1 

H = EffiW. v = i, 

i=l j=0 

for e = 0, 1; (2.7) 

■y n k—l 

V k = (x n+1 + B)V k - 1 + -Y^iT, V k-i- j el V = 0, 

1 i=l j=l 

for e = A{u- x n+1 + B). (2.8) 
The explicit formulae for the integrals Hi have the form 

n I M 2 N 

H t = -E^^ + ^ + ^EV fc ef- fc , for £ = 0,1, (2.9) 
— p ■ — p ■ — 

i=i c * C J fc=o 

AT-l 

# 4 = 2x 2 Y {- l ) 3 ^4 Vn-j + 4p n+ i Pi a:< - P- (e* + 4x n+1 - 45) 
i=i 
M 2 

+ E^^' for e = 4(u-x n+1 + 5), (2.10) 
where Mij = XiPj — XjPi. The Hamiltonians H are given by 

n n 

H = YH l = e-YP* + V N+ i, for £ = 0,1, (2.11) 

1=1 i=l 
I n+1 

H = ~E^ + V ^ for e = 4(u-x n+1 + B). (2.12) 

1 8=1 

Proof is straightforward and based on direct computations. 

We remark, that the above recurrence formulae for the potentials can be written in 
differential form. In particular, for the paraboloidal coordinates we have 

9Vat = l dVjv-i ^ _ A dVx~\ ^ i = i n (2.13) 

OXi 2 OX n+ i OXi 

dV N 1 " dV N ^ . ^QVn-i 



Viv-i + ■= E^^P + K+i - • (2.14) 



OT n+ i 2 dxi ax n+ i 
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Notice that the case of e = is connected with the ellipsoidal coordinates on a sphere 
and two other cases e = 1, 4(ii — x n+ \ + B) describe the ellipsoidal and paraboloidal 
coordinates in the Euclidean space, respectively (see Section 3.2 and [11], for more 
details). Recall that we study now the motion of a particle on these manifolds under 
the external field with the potential V that could be any linear combination of the 
homogeneous ones 

Now we are ready to describe the linear algebra for the L-operator (|2.1|). 



Theorem 2 Let the L-matrix be of the form j\2. J\ ) and satisfy the conditions of the 
Theorem 1, then the following algebra is valid for its entries 



{B(u),B(v)} 
{C N (u),C N (v)} 

{B(u),A(v)} 
{C N (u),A(v)} 
{B(u),C N (v)} 



{A(u),A(v)} = 0, 
-ia N (u,v) (A(u)-A(v)), 

2 -(B(u)-B(v)), 



u — V 
2 

u — V 
4 



u — V 

where the function ctjs[(u,v) has the form 

Qn(u) - Qn(v) 



-(C N (u) - C N (v)) -2a N (u,v)B(v) 
(A(u)-A(v)), 



cxn(u, v) 



Qn{u) 



u — V 



N 



k=l 
k 



k k 
U — V 



u — V 



(2.15) 
(2.16) 

(2.17) 
(2.18) 
(2.19) 

(2.20) 



^2QkU k , Qk — ^2 VmVk-m- 
fc=0 m=0 



Proof is based on the recurrence relations 

We remark that for the paraboloidal coordinates the following formula is valid 



Q(u) 

therefore in this case we have 



7V-2 



Q(u) 



1 ^ Qy^^ 

2 k=0 dx n+1 
1 dC{u) 



(2.21) 



(2.22) 



4 dx n+1 ' 

The algebra (|2.15| )- (|2.19| ) can be rewritten in the matrix form as linear r- matrix 
algebra 

{Lx{u),L2(v)} = [r(u - v), L x {u) + L 2 (v)\ + [s N (u } v), L^u) - L 2 (v)\, (2.23) 

/, L 2 {v) = I (g) L(v); the matrices r{u — v ) and 



using 4x4 notations Li(u) = L(u) 
Sn(u,v) are given by 



r(u — v) 



u — V 



p, p 



s N (u, v ) = 2 a N (u, v) cr_ ® cr_, 



^ 1 ^ 
10 
10 

V o o o i J 

' 
, 1 o 



(2.24) 



The algebra ( 2.15|) -( 2.19|) or (|2.23|) - (|2.24 ) contains all the information about the system 
under consideration. From it there follows the involutivity of the integrals of motion. 
Indeed, the determinant d(u) = detL(tt) is the generating function for the integrals of 
motion and it is simply to show that 

{d{u),d{v)} = 0. (2.25) 

In particular, the integrals Hi ( |2.9|) , (|2.10|) are the residues of the function d{u): 

Hi = res| u=ei d(u), i = l,...,n. 

The Hamiltonians H (|2.11|) ,( 2.12 ) appear to be a residue at infinity. Let us rewrite the 
relation ( [2.231 ) in the form 



with dij — Tij -\- Sij } dji — Sij 



{L 1 (u),L 2 {v)} = [di 2 (u,v),Li(u)} - [d 2 i{u,v)L 2 {v)}, 
Tij at i < j. 



(2.26) 



Theorem 3 The following equations (dynamical Yang-Baxter equations) are valid for 
the algebra (\2.26j) 

[di 2 (u, v), dx 3 (u, w)] + [d 12 (u, v), d 23 (v, w)] + [d 32 (w, v), d 13 (u, w)] + 
+{L 2 (v), d i3 (u, w)}- {L 3 (w), d i2 (u, v )} + 

+[c(u, v, w), L 2 (v) - L 3 {w)\ = 0, (2.27) 

where c(u,v,w) is some matrix dependent on dynamical variables. Other two equations 
are obtained from ($.2% ) by cyclic permutations. 

Proof Let us write the Jacobi identity as 



{Li(u), {L 2 (v) , L 3 (w)}} + {L 3 (w), {Li(u), L 2 (v)}} 
+{L 2 (v),{L 3 (w),L 1 (u)}} = 



(2.28) 



with Li(u) = L(u)®I®I, L 2 (v) = I®L(v)®I, L 3 (w) = I®I®L(w). The extended 
form of the ( [2.28D reads [13|, 



[Li(u), [di 2 (u,v),di 3 (u,w)] + [d 12 (u,v),d 23 (v,w)] 
+[d 32 (w,v),d 13 (u,w))) + 

+ [Li(u), {L 2 (v),d 13 (u,w)} - {L 3 (w),d 12 (u,v)}] - 
+ cyclic permutations = 0. 



(2.29) 



Further we restrict ourselves proving the ( |2.27|) only in the paraboloidal case (other 
cases can be handled in a similar way). To complete the derivation of (|2 . 2 7| ) we shall 
prove the following equality for all the members of the hierarchy 



{L 2 (v),s 13 (u,w)} 



{L 3 (w),s 12 (u,v)} = 2f3 N (u,v,w){P 23 ,s 13 + s 12 ] 
8Pn(u, v, w) 



3x n 



+1 



-[s,L 2 (v) - L 3 (w)] 



(2.30) 
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(with cyclic permutations). In ( |2.30|) the matrix s — er_ (g) cr_ ® (X_ and 

Qn(u)(v -w) + Qn(v)(w -u) + Q n (w)(u - v) 



/3 N (u,v,w) 



(u — v ) (v — w) (w — u) 
In the extended form ( |2.30[ ) can be rewritten as 
{Q(u),Q(v)} 
{A(u),Q(v)} 



{B(u),Q(v)} = 0, 

1 da]y(u, v) 



2 dx n+ i 



B(u), 



{A(w),a N (u,v)} 



(a N (w, u) - a N (w, v)) 

u — v 

\B{w) (da>N{w,v) dajy(w,u) 
2 u — v \ dx n+ i dx 



n+l 



{Q(u),C(v)} 



{C(u), Q(v)} = \ A(u) - A{v)). 

OX n+ i 



(2.31) 

(2.32) 
(2.33) 

(2.34) 
(2.35) 



The equality (|2.32|) is trivial and equation ( p.33|) is derived by differentiating ( p,18|) . 
Equation ( |2.34j ) follows from the definition of Q(u) and ( ^.33 ). To prove ( 2.35 ) we write 
it using the explicit form of Cn{u) and A(u) as follows 



E 



Pi 



\ (v - ei)(u - e { ) 
1 ^--^ i>j.i i 



(u 



dajsiiw, u) 



iv - e, ; 



dx,i 



da N (w,vY 
dxi 



2 ( u ~ e i)( v - e i 
Using the identity 



dx 



k h 

w — u 



u- 



w — u 



h h 

w — V 



w — V 



n+l 



w k+1 — u k+1 



W — U 



dx n 



(2.36) 



+i 



w 



k+l 



,fc+l 



W — V 



and the recurrence relation (|2.13|) , we find that the equality ( p.36| ) is valid. Therefore 
the equations ( |2.27| ) follow with the matrix c(u, v, w) = dj3(u, v, w)/dx n+ i a_ (g) o"_ (g) a_ . 
The proof is completed. 

We remark that validity of the equations (|2.27|) with an arbitrary matrix c(u, v, w) 
is sufficient for the validity of ( j2.28|) and, therefore, ( |2.27|) can be interpreted as some 
dynamical classical Yang-Baxter equation, i.e. the associativity condition for the linear 
r-matrix algebra. These equations have an extra term [c, Li — Lj] in comparison with 
the extended Yang-Baxter equations in [|l]| . 

We would like to emphasize that all statements of the Section can be generalized to 
the following form of the potential term V/v(w) in 



VmnM 



N 

E 

k=-M 



fkV k u 



N-k 



f k eC 



that corresponds to the linear combinations of homogeneous terms 14 as potential V 
and also includes the negative degrees to separable potential. See the end of Section 
3.1 for more details. 
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3 Consequences of the r-matrix representation 
3.1 Lax representation 

Following the article || we can consider the Poisson structure ( |2.26 ) for the powers of 
the L-operator 



{(L 1 («))^(L 2 (^)) ^ }=[^ ) ( W ,^;),L 1 («)]-[4f ) («^)^2(^)], (3.1) 



with 



4 r } (u,v)= (3.2) 



»(*.«), 

Hj 

EE(ilW)^ 1 (i2(t ) )) M - 1 rf lJ («^)(i 1 (u)) !, (i2(«)) 9 . 
p=0 g=0 



As an immediate consequence of ( |3.1|) - (|3.2|) we obtain that the conserved quantities H, 
Hi are in involution. Indeed, we have 

{Tr(L 1 ( M )) 2 ,Tr(L 2 (t;)) 2 } = Tr{(L 1 ( M )) 2 ,(L 2 (t;)) 2 }, (3.3) 



and after applying the equality ( |3.1| ) at k — I — 2 to this equation and taking the 
trace, we obtain the desired involutivity. Further, let us define differentiation by time 
as follows: 

L(u) = jL{u) = Tr 2 [L 1 (u) ) (L 2 (u)) 2 }, (3.4) 

where the trace is taken over the second space. Applying the equation ( |3.1| ) at k = 
1, / = 2 to Q3.4|) , we obtain the Lax representation in the form L(u) = [M(u), L(u)] 
with the matrix M(u) given as 

M(u) = 2\im r Ti 2 L l (v)(r(u - v) - s(u,v)). (3.5) 

After the calculation in which we take into account the asymptotic behaviour of the 
L-operator fl2.1|) , we obtain the following explicit Lax representation: 



L(u) = [M(u),L{u)\, 

where Qn(u) was defined by the equations ( |2.20| ). Lax representations for the higher 
flows can be obtained in a similar way. 
It follows from the equation (|3.6p that 



A(u) = -\b(u), C n {u) = - l -B{u) - B(u)Q N {u), 



2 v ' y ' 2 

so our L-matrix can be given in the form 



-("§?(«) fi H ^ (3 7) 
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The equations of motion, which follow from ( |3.6| ) with the L- matrix from (|3.7|) , have 
the form 

B 1 [Q N ]-B(u) = 0, (3.8) 

where 

B x [Q N ] = - A d* + ^{d,Q N }, d=j t , 

with curly brackets standing for the anticommutator. Operator B\ is the Hamiltonian 
operator of the first Hamiltonian structure for the coupled KdV equation ffl, |2J. Equa- 
tion ( |3.8j ), considered as one for the unknown function B(u), was solved in the three 
cases 

n x 2 



B(u)=e-J2 — — > e = 0,l,4(u-x n+1 +B), 
in |L| and 0. General solution of this equation as one for the Q(u) has the form 

JV 

Qmn(u)= fkQkU k , f k eC, (3.9) 

k=-M 

where the coefficients Qk are defined from the generating function Q(u) 

+oo 

Q(u) = B' 2 (u) = QkU*- (3-10) 

k=— oo 



Recall that we can write the element Cn(u) of the L- matrix ( [2.1|) in two different forms 
(using Q or V functions) 

1 •■ n v 2 

C N {u) = --B(u) - B(u)Q N (u) = - - V(u) 

where function V(u) = Y, k =oVkU N ~ k was defined in ( |2.4| ). The general form of the 
function V{u) is as follows: 

N 

Vmn(u)= J2 fkV k u k , f k eC, (3.11) 

k=-M 

where coefficients V k are defined by the generating function V(u) 

+oo 

V(u) = B- l (u) = V kU k . (3.12) 

k=— oo 

Potentials V k are connected with coefficients Qk- Indeed, using generating functions 
( CT ) and (|3~T2D , we have 

(+oo \ / +oo \ +oo k 

fc=— oo y \fc=— oo y fc=— oo i=o 



and, therefore, = J2j=o^k^k~j- Thus we have recovered the formula ( |2.20| ) for the 



•.k 

s-matrix. 
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3.2 Separation of Variables 



Let K denote the number of degrees of freedom: K = n — 1 for ellipsoidal coordinates 
on a sphere, K = n for ellipsoidal coordinates in the Eucludean space, and K = 
n + l for paraboloidal coordinates in the Euclidean space. The separation of variables 
(c.f.[]T2"|, [T^ D is understood in the context of the given hierarchy of Hamiltonian systems 
as the construction of K pairs of canonical variables 7r», fa, i — 1, . . . , K, 

{l^h A*fc} = {^i, TTfc} = 0, {iri, jj k } = 6ik, (3.13) 
and K functions $j such that 



(i) 

N 5 



HP) 



0, J = 1,2, 



(3.14) 



where are the integrals of motion in involution. The equations (|3.14 ) are the 
separation equations . The integrable systems considered admit the Lax representation 
in the form of 2 x 2 matrices ( |3.6|) and we will introduce the separation variables 7Tj, fa 
as 

B(ni) = 0, TCi = A(iMi), i = l,...,K. (3.15) 

Below we write explicitly these formulae for our systems. The set of zeros fa,j = 1, ... K 
of the function B(u) defines the spherical (e = 0), general ellipsoidal (e = 1) and general 
paraboloidal (e = 4(w — x n+1 + B)) coordinates given by the formulae B, p, |T! 



/;/ 



where c = V] 



x 



k- 



1,- 

for e 



0; 



fc=i 



;r; 



^n+1 



-4 



- e fc J 
for e — 1; 

n n+l 
- ^ e i + -6 + 53 ^ 



(3.16) 



(3.17) 



i=i i=i 
-4- 



for £ = 4(u — X n+1 + S). 



(3.18) 



Theorem 4 T/ie coordinates fa, 7Tj ^iu en fry ( \3.13j ) are canonically conjugated. 
Proof Let us list the commutation relations between B(v) and A(u), 

{B(u),B(v)} = {A(u),A(v)} = 0, (3.19) 

{A(u), B(v)} = —(B{u)-B{v)). (3.20) 
v — it 
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The equalities = follow from ( ft.lSQ . To derive the equality {fi^iVj} = —5ij 

we substitute u = jij in ( |3.20| ), obtaining thus 



-B(v) 



which together with the equation 

= {lTj,B(Hi)} = {TTj,B(v)} \ v=fli +B'(fl i ){lT j ,fX i } 

gives 



{TTjj/ii} 



1 



B'(jn 



-{Kj,B(v)} \ V=IH = Sij. 



Equalities {7r i; ttj} = can be verified in the similar way: 

{7ii,7ij} = {A(m),A(nj)} 

= {^O;M0)} l«= w +A'(fi l ){ f i i ,A(fi j )} 
= A'(fj,j){A(fii),fj,j)} + A'iinXfH, Afjij)} = 0. 

The separation equations have the form 

vif = d(fii), 

where the function d(u) is the determinant of the L-operator (|2.25 ). 



(3.21] 



3.3 Quantization 



The separation of variables has a direct quantum counterpart [11, |2(J. To pass to 
quantum mechanics we change the variables /ij to operators and the Poisson brackets 
( |3.13j ) to the commutators 



\}ij,Hk] = [7Tj,7r fe ] = 0, [nj,fi k ] = -iS jk . 



(3.22) 



Suppose that the common spectrum of \Xi is simple and the momenta 7Tj are realized 



as the derivatives nj = The separation equations (|3.21|) become the operator 



equations, where the noncommuting operators are assumed to be ordered precisely in 
41 * u ci i> <i->* ; -- - - H$, . . .,Hj^\ Let ^(/ii, . . . , Hk) be 



the order as those listed in ( p. 14 ), that is Ki,fii,H 



a common eigenfunction of the quantum integrals of motion: 

fl$tt = Aitt, i = l,...,K. (3.23) 
Then the operator separation equations lead to the set of differential equations 

,K, (3.24) 
which allows the separation of variables 

K 

(3.25) 



*(/*!,..., H K ) = Y[lpj(fMj) 

3=1 
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The original multidimensional spectral problem is therefore reduced to the set of one- 
dimensional multiparametric spectral problems which have the following form in the 
context of the problems under consideration 

for e = 0,l, (3.26) 

/ d 2 n A \ 

— + l6u N - 2 (u + B) 2 + 8A n+1 + J2 —T~ ^ Al ' • • • ' A ™+!) = °' 

for e = A(u-x n+1 + B), (3.27) 

with the spectral parameters Ai, . . . , X n +i- The problems ( p.26|) , (|3.27|) must be solved 
on the different intervals — "permitted zones" — for the variable u. 



4 Conclusion 

We remark that all systems considered yield the algebra which has general properties 
being independent on the type of the system. Therefore it would be interesting to 
consider its Lie-algebraic origin within the general approach to the classical r-matrices 
lj. 

There exists an interesting link of the algebra studied here with the restricted flow 
formalism for the stationary flows of the coupled KdV (cKdV) equations Q. The 
Lax pairs which have been derived in the paper from the algebraic point of view were 
recently found in by consideration of the bi-Hamiltonian structure of cKdV. 

It seems to be interesting to study the same questions for the generalized hierarchy 
of differential operators of Gelfand-Dickey for which the corresponding L-operators have 
to be the n x n matrices. 
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